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Abstract 

We revisit the role of attractions in liquids and apply these concepts to colloidal suspensions. 
Two means are used to investigate the structure; the pair correlation function and a recently 
developed topological method. The latter identifies structures topologically equivalent to ground 
state clusters formed by isolated groups of 5 < m < 13 particles, which are specific to the system 
under consideration. Our topological methodology shows that, in the case of Lennard-Jones, the 
addition of attractions increases the system's ability to form larger (m > 8) clusters, although pair- 
correlation functions are almost identical. Conversely, in the case of short-ranged attractions, pair 
correlation functions show a significant response to adding attraction, while the liquid structure 
exhibits a strong decrease in clustering upon adding attractions. Finally, a compressed, weakly 
interacting system shows a similar pair structure and topology. 

PACS numbers: 82.70.Dd; 82.70.Gg; 64.75. +g; 64.60.My 
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I. INTRODUCTION 



Among the cornerstones of our understanding of the structure of bulk simple liquids is 
that it is dominated by the repulsive core. This leads to the idea that hard spheres form a 
suitable basic model of the liquid state. The liquid pair structure may then be accurately 
calculated using, for example the Percus-Yevick closure to the Ornstein-Zernike equation for 
hard spheres, and treating the remainder of the interaction as a perturbation^-. 

Although in principle colloidal dispersions are rather complex multicomponent systems, 
the spatial and dynamic asymmetry between the colloidal particles (10 nm-1 /xm) and smaller 
molecular and ionic species has enabled the development of schemes where the smaller 
components are formally integrated out 4 . This leads to a one-component picture, where 
only the effective colloid-colloid interactions need be considered. The behaviour in the 
original complex system may then be faithfully reproduced by appealing to liquid state 
theory^ and computer simulation 6 . Since the shape of the particles is typically spherical, 
and the effective colloid-colloid interactions may be tuned, it is often possible to use models 
of simple liquids to accurately describe colloidal dispersions. 

In colloidal systems, due to the mesoscopic length- and longer time-scales, one may 
also determine the structure in real space in 2D and 3D at the single particle level using 
optical microscopy- 1 - and optical tweezers^. This may be done with sufficient precision that 
interaction potentials can be accurately determined both for purely repulsive systems^ 1 ^ 
and for systems with attractive interactions^ 1 -. 

It has been conjectured as far back as the 1950s that the structures formed by clusters 
of small groups of particles in isolation might be prevalent in liquids 1 ^ 1 ^. More recently it 
has been demonstrated that for spherically symmetric attractive interactions, the structure 
of clusters of size m > 7 particles depends upon the range of the potential, as shown in 
Fig. [I] 14 . This brings a natural question: if the structures defined by these clusters are 
indeed prevalent in liquids, and they depend upon the range of the interaction, then might 
liquids with differing interaction ranges exhibit differing cluster populations? Moreover, 
while removing the attractive component of Lennard- Jones [Fig. [2(b)] has little effect on 
the pair structure 1 -^, what is the effect on any cluster population? 

We have recently developed a novel means to identify structure in simple liquids. In 
isolation, small groups of particles form clusters with well-defined topologies. These have 
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Figure 1: (color online) Clusters found in bulk systems using the topological cluster classification. 
For m < 7, where m is the number of particles in a cluster, all studied ranges of the 
Morse potential equation (J3j) form clusters of identical topology. In the case of larger m the cluster 
topology depends on the interaction range. Here we follow the nomenclature of Doye et. alM- where 
A corresponds to long-ranged potentials and B.... to minimum energy clusters of shorter-ranged 
potentials. 

been identified for the Lennard- Jones potential^ and for the Morse potential, which has 
a variable ranged. We identify clusters relevant to the Lennard- Jones and Morse Poten- 
tials in bulk liquids, with a method we term the Topological Cluster Classification (TCC)^. 
Here we use this scheme as a highly sensitive probe of the liquid structure. It is our inten- 
tion to use the TCC to investigate possible differences in structure between the Lennard- 
Jones liquid and that resulting from the repulsive part of the Lennard- Jones interaction, the 
Weeks- Chandler- Andersen (WCA) potential^. Although we have argued that some colloidal 
liquids are well described by a short-ranged Morse potential, the structure of clusters 
of adhesive hard spheres has very recently been shown to exhibit some degeneracy, with 
multiple cluster topologies having the same number of bonds in the limit of short-ranged 
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Figure 2: (color online) Interaction potentials used, (a) Long-ranged potentials: Morse (dark green) 
and truncated Morse (bright green) with range parameter po = 4.0. (b) Lennard-Jones (red) and 
WCA (pink), (c) Short-ranged potentials: Morse (blue) and truncated Morse (turquoise) with 
range parameter po = 25.0. Dashed cyan line in (c) denotes the hard sphere interaction, oeff 
denotes the effective hard sphere diameter as defined in equation ([7]) and listed in tabled 

attractions^ 1 ^. However, minimising the second moment (or radius of gyration) of clusters 
of hard spheres 2 ^" 2 ^ shows a strong similarity with the short-ranged Morse system 14 . Our 
TCC methodology has some similarities to the common neighbour analysis introduced by 
Andersen22>22, however here we focus on clusters rather than bonds. 

Since the tunability of colloidal systems allows a wide range of potentials to be realised, 
including long-ranged interactions relevant to metals^^ 4 - 1 ^, we also consider long-ranged 
(Morse) potentials, in addition to the Lennard-Jones interaction and short-ranged Morse 
potential, along with their purely repulsive counterparts. We further compare with hard 
spheres. In these systems, we study the groups of particles topologically equivalent to 
ground state clusters found in isolation. 

The canonical model of colloid-polymer mixtures of Asakura and Oosawa, assumes hard 
sphere colloid-colloid and colloid-polymer interactions, while the polymer-polymer interac- 
tion is ideal 26 . A one-component description^ 1 ^, accurate for small polymer-colloid size 
ratios 2 ^ leads to a hard core with a short-range attraction. We have recently shown that, 
for the parameters we shall consider here, the continuous Morse potential provides a reason- 
ably accurate description of this system 2 ^. Meanwhile, longer interaction ranges correspond 
to metals^ and purely repulsive long-ranged interactions are relevant to soft matter sys- 
tems such as charged colloids, star polymers 3 -, star polyelectrolytes^l and colloidal microgel 
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particles. 



This paper is organised as follows. In section [TT] we describe the simulation methodology 
and our approach for comparing different interaction potentials, our results are presented 
in section II III and we conclude with a discussion in section IIVI Our main results can be 
summarised as follows. Although Lennard- Jones shows very little change in the radial dis- 
tribution function g(r) upon adding attractions, the topology is significantly altered: adding 
attractions promotes the formation of larger clusters. Conversely, short-ranged systems show 
the opposite behaviour: adding attractions strongly decreases clustering, while the first peak 
of g(r) shows some increase upon adding attractions. 

II. SIMULATIONS AND INTERACTION POTENTIALS 

We use standard Monte-Carlo (MC) simulation in the NVT ensemble^ with N = 2048 
particles. Each simulation run is equilibrated for 2 x 10 7 MC moves and run for up to a 
further 10 s moves. In all cases, we confirmed that the system was in equilibrium on the 
simulation timescale by monitoring the potential energy. 

A. Interaction potentials 

We seek to compare systems with different interactions, under similar conditions. Weeks, 
Chandler and Andersen^ provided a protocol by which the Lennard- Jones potential could 
be compared with the so-called WCA potential (Lennard- Jones without attractions). The 
pair interaction u(r) is separated into two parts: 



where r is the separation between particles, Uq(t) is the reference (repulsive) interaction and 
w(r) is the perturbative attraction. In the Lennard- Jones case, 



where (3 = l/Zc^T where kg is Boltzmann's constant and T is temperature. Here Elj = 1/T 
is the well depth. WCA thus define the reference potential as 



u{r) = u (r) + w{r) 




(1) 
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(-) -(-) +(3e LJ forr<2 1 /V, 
Pu WCA (r)={ ' UJJ " (2) 

for r > 2 1 / 6 a. 

2nd order perturbation theories 1 allow accurate prediction of the pair structure. However, 
here we are interested in a particle based analysis that probes the structure at a level 
beyond the two body distribution function and restrict ourselves to the interactions given 
by equations (JTJ and (j2J). 

In the case of the longer and shorter ranged interactions, we use the Morse potential 
which reads 



f3u M (r) = Pe M e p ^- r) (e P0 ^ - 2) (3) 

where po is a range parameter and (3em is the potential well depth. We set po — 25.0 
to simulate a system with short-ranged attractions similar to a colloid-polymer mixture 
and po — 4.0 as an example of a longer-ranged system. Following the WCA approach, we 
introduce a repulsive (truncated) Morse potential, which is also truncated at the minimum 
and is defined as follows. The potentials we use are plotted in Fig. [2j 

f f3e M e p0 (<T " r) (e p0 (<T - r) - 2) + (3e M for r < a, 
(3u T M(r) = < (4) 

I for r > a, 

This truncated Morse potential is thus similar to hard spheres for p = 25.0 [Fig. [2(c)]. 
The repulsive systems have well-defined truncations, 2 1 / 6 cr for WCA and a for the Morse 
potential. In the case of the attractive systems, we truncate and shift both Lennard- Jones 
and Morse (po = 25.0) at 2.5a and the long-ranged Morse (po = 4.0) at 4. Oct. 



B. Comparing different systems 

We have outlined a means by which we can compare systems with and without attraction, 
by removing the attractive part of the interaction. In order to match state points between 
systems with differing interaction ranges, we use the extended law of corresponding states 
introduced by Noro and Frenkel 32 . Specifically, this requires two systems to have identical 
reduced second virial coefficients B\ where 
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B *2 = xJf- (5) 



S^CTeFF 



where ueff is the effective hard sphere diameter and the second virial coefficient 



oo 



B 2 = 2n J drr 2 [1 - exp (-/3u{r))] . (6) 

Q 

The effective hard sphere diameter is defined as 

oo 

oeff = J dr[l- exp (-(3u RE p(r))} (7) 
o 

where the repulsive part of the potential urep is described above in section III Al Thus we 
compare different interactions by equating B\ and ueff- The latter condition leads to a 
constraint on number density 

Nn(JEFF /m 

PEFF = — Qy- ( 8 ) 

where V is the volume of the simulation box. We fix Peff to a value equivalent to the 
Lennard- Jones triple point (pu = 0.85) throughout. In the case of hard spheres, this value 
is pus ~ 0.9310 or (pus = t^Phs/Q ~ 0.4875 where is the packing fraction. Details of state 
points investigated are given in table HI 



C. The Topological Cluster Classification 

To analyse the structure, we first identify the bond network using a modified Voronoi 
construction with a maximum bond length r c = 1.8a and four-membered ring parameter 
f c = 0.82 (TCC) paper—. Having identified the bond network, we use the TCC to determine 
the nature of the cluster. This analysis identifies all the shortest path three, four and five 
membered rings in the bond network. We use the TCC to find clusters which are global 
energy minima of the Lennard- Jones and Morse potentials. We identify all topologically 
distinct Morse clusters, of which the Lennard- Jones clusters form a subset (the Lennard- 
Jones and Morse interactions are similar in the case that the range parameter po = 6.0 and 
the topology of the ground state clusters is identical). In addition we identify the FCC 
and HCP thirteen particle structures in terms of a central particle and its twelve nearest 
neighbours. We illustrate these clusters in Fig. [TJ In the case of the Morse potential, for 
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m > 7 there is more than one cluster which forms the ground state, depending on the 
range of the interaction^. We therefore consider ground state clusters for each system and, 
separately, calculate all topologically distinct Morse clusters for m < 14. For more details 
se©i£. 

To compare the various fluids we study here, we proceed as follows. Comparing systems 
with and without attractions, we consider the ground state clusters (of the attractive sys- 
tem). If a particle is a member of more than one cluster, it is taken to 'belong' to the larger 
cluster. Thus, the total cluster population < iV the total number of particles. However, 
when we seek to compare different potentials, we need to account for the fact that these 
may have different ground state clusters. If the particle is part of two clusters which are 
different in size, we choose to count it as the larger cluster, but if the particle is part of two 
clusters of the same size, it is counted as the cluster corresponding to the shorter-ranged 
interaction. In this case, the total number of particles counted as belonging to a cluster can 
exceed the number of particles in the simulation. 

D. Systems studied 

The different systems considered are listed in Table [H In addition to the state point 
(e,p), we list the reduced second virial coefficient B\ and effective hard sphere diameter 
ueff- Some comments upon the use of clusters in the case of repulsive systems are in order. 
Clearly, isolated clusters require cohesive forces, however here we compare the WCA and 
the truncated and shifted Morse potential with their cohesive counterparts and we assume 
it is appropriate to consider the same clusters. Given the similarity of the truncated Morse 
potential to hard spheres [Fig. [2(c)], it is instructive to include these also. 

We are motivated to consider the Lennard- Jones triple point, as we expect clusters to be 
more prevalent at lower temperature. However, mapping the short-ranged Morse potential 
(po = 25.0) to the Lennard- Jones triple point leads to a system unstable to crystallisation. 
We found that at higher temperature the Morse system was stable (on the timescales of 
these simulations) against crystallisation for (3em = 2.0. Thus we compare high-temperature 
Lennard- Jones (T=2.284) with Morse po = 25.0 and triple point Lennard- Jones with the 
long-ranged Morse po = 4.0. 
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Table I: State points studied. LJ high temp, and triple correspond to the two temperatures at 
which Lennard-Jones and WCA simulations were carried out. Trunc. Morse denotes the truncated 
Morse interaction (equation [3]) . 



System 


B* 


08 


p °EFF 


LJ high temp. 


-0.2325 0.4447 0.9776 0.9839 


WCA high temp. 


2.013 


0.4447 0.9776 0.9839 


Morse p = 25.0 


-0.2325 


2.0 


0.9837 0.9818 


Trunc. Morse po = 25.0 


0.9818 


2.0 


0.9837 0.9818 


Hard Spheres 


1.0 


N/A 


0.9310 1.0 


LJ triple 


-3.742 


1.471 


0.85 1.0308 


WCA triple 


2.307 


1.471 


0.85 1.0308 


Morse p = 4.0 


-3.742 


0.9109 


1.548 0.8441 


Truncated Morse po = 4.0 


1.015 


0.9109 


1.548 0.8441 



III. RESULTS AND DISCUSSION 

A. The Lennard-Jones Triple Point: Long-ranged interactions 

We take as our starting point for the analysis of these data the result that for dense 
liquids, the WCA potential readily captures the pair structure of the Lennard-Jones liquid^. 
The radial distribution functions g(r) for the Lennard-Jones liquid at the triple point and 
the corresponding WCA system are plotted in Figure [3(b) . The effectiveness of WCA in 
describing the pair structure is clear. The same observation holds for the longer-ranged 
(p = 4.0) Morse and truncated Morse systems shown in Figure E(a). 

Turning to the TCC analysis, in the WCA-Lennard- Jones system we see a somewhat 
different story, as shown in Fig. 13(b). Note the logarithmic scale in this plot: cluster 
populations vary over three orders of magnitude, clear relative differences are seen between 
Lennard-Jones and WCA. Unlike the g{r) which are very similar, there is a clear trend 
in the cluster populations N c /N. Larger clusters are more prevalent in the Lennard-Jones 
system, compared to the WCA. The differences are emphasized in Fig. 0(b) which plots the 
ratio of the cluster populations in Fig. 11(b). One might argue that smaller clusters may 
readily be formed simply by compressing spheres. The addition of attractions promotes the 

9 




Figure 3: (color online) Pair-correlation functions, (a) Long-ranged potentials: Morse po = 4.0 
with (dark green, dashed) and without (bright green) attractions. Here /3em = 0.9109. (b) 
Lennard- Jones (red, dashed) and WCA (pink) for a well depth of = 1-471 (the triple point). 

formation of larger clusters which require more organisation and co-operativity. We remark 
that the difference in structure revealed by the TCC is rather significant, given that the 
radial distribution functions are so similar. For example, there is a twofold difference in 
the triangular bipyramid 5A, one of the most popular clusters. As for the 13A icosahedron, 
its population is quadrupled by adding attractions. Note also that in these equilibrium 
liquids, we find a small but measurable number of particles with local crystalline topology, 
even though there is no sign of splitting in the second peak of g(r) (Fig. [3]), which is often 
taken to be a sign that the liquid is close to crystallising—. Of the smaller clusters, the 7A 
pentagonal bipyramid is found in limited quantities. However, it is found also as part of all 
the larger clusters which form minima for the Lennard- Jones except HCP and FCC so our 
counting methodology counts some 7A (D 5 h) particles as members of larger clusters. 

The longer-ranged Morse (po = 4.0) system on the other hand shows very little difference 
upon adding attractions. Due to its softness, the long-ranged Morse system has a rather 
small value of ueff (table H]). Thus, matching Lennard- Jones requires a higher density, 
which leads to some overlap of the particles. For example the mean inter-particle spacing 
dm = p" 1 / 3 ~ 0.8644a. The Morse potential has its minimum located at a. However, such is 
the long ranged nature of Morse po — 4.0 that in fact, 0.8644a remains within the attractive 
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Figure 4: (color online) Population of particles in a given cluster. iV c is the number of particles 
in a given cluster, N the total number of particles sampled. Here we consider only ground state 
clusters for each system, (a) Morse (po = 4.0) (dark green) truncated Morse (bright green), (b) 
Lennard- Jones at the triple point (red) and corresponding WCA (pink). Note the semi-log scale. 

well, although there is some compression. Thus both with and without attractions, the 
system is compressed, which may dominate the local structure. Furthermore, the value of 
@Em = 0.9109 (table HJ indicates that the interactions here are rather weak. As we shall see 
below, weaker interactions can lead to topologically similar structures. Note the relatively 
high abundance of the pentagonal bipyramid 7A, due to fewer reclassifications as higher- 
order clusters. Plotting the ratio of the cluster populations [Fig. Ufa)] further emphasizes 
the similarity of the long-ranged Morse interaction with and without attractions: the ratio 
lies close to unity in all cases, except the 13A icosahedron where the statistics are insufficient 
for a reliable comparison. 

We now plot the population of all identified clusters. This enables us to directly compare 
the populations of the long-ranged Morse and Lennard- Jones systems. We see that there 
is no strong preference for ground states [shaded in Fig. (jB])]. In fact a number of ground 
states are less populated than other clusters of the same size. Comparing the different 
systems, the general trend is of LJ/WCA tending to form larger clusters than the long- 
ranged Morse, which is consistent with the idea that the long-ranged Morse is a weakly 
interacting, compressed fluid. Recall that, for example in the case of 11-membered clusters, 
we count a given particle as a member of an 11F (C% v ) if it is a member of more than one 
m — 11 clusters of which one is an 11F. While this may innate the populations of such 
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5A 6A 7A 8A 9A 10B 11C 12B 13A HCPFCC 5A 6A 7A 8B 9B 10B 11C 12B 13A HCPFCC 

Cluster Type Cluster Type 

Figure 5: (color online) Ratio of cluster populations in systems mapped to the Lennard-Jones triple 
point, (a) Morse and truncated Morse (po = 4.0) and (b) Lennard-Jones and WCA. These plot 
the same data as Fig. d] expressed to emphasize the difference between the systems. 

clusters, we argue that systems with differing ground states are compared in an unbiased 
way. One result of considering all Morse clusters is that Lennard-Jones has by far the largest 
number of 13A icosahedra, although the population of the 13B decahedral cluster is larger. 

B. High- Temperature systems : Short-ranged interactions 

For shorter-ranged interactions relevant to colloid-polymer mixtures^, to avoid crystalli- 
sation we used an attractive well depth of /3em = 2.0, which corresponds via equation (J6j) 
to a Lennard-Jones well depth of (3elj ~ 0.4447. Pair correlation functions are shown in 
Fig. [71 Again, the WCA and Lennard-Jones [Fig. [TJa)] show a similar behaviour. In the 
case of the shorter-ranged potentials [Fig. Wip) and inset], we see a strong increase in the 
first peak. The short-ranged Morse po — 25.0 system shows some splitting of the second 
peak. We carefully checked that no crystallisation was found during these simulation runs, 
which is supported by the TCC analysis which shows a much reduced population of par- 
ticles in a locally crystalline environment upon adding attractions (Fig. |8]). However, we 
are unaware of an equilibrium phase diagram for the Morse po = 25.0 system, so we cannot 
exclude the possibility that the system is metastable to crystallisation, a point to which we 
return below. Furthermore, the first peak in g(r) is rather higher in the case of the Morse 
interaction, compared to the purely repulsive truncated Morse and hard-sphere interactions. 
We note also that there is little difference between the truncated Morse and hard sphere 
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Cluster Type 

Figure 6: (color online) Population of particles in a given cluster at parameters mapped to the 
Lennard-Jones triple point. N c is the number of particles in a given cluster, N the total number of 
particles sampled. Here we consider ground state clusters for all ranges of the Morse potential^. 
Colours are Lennard-Jones (red), corresponding WCA (pink), Morse (po = 4.0) (bright green) and 
truncated Morse (dark green). Those clusters which are ground states are labelled as 'both' when 
both potentials share the same ground state, and 'LJ' and 'M' corresponding to the Lennard-Jones 
and Morse cases accordingly. Note the semi-log scale. 

pair correlation functions. This suggests that the truncated Morse po — 25.0 may provide 
a useful continuous approximation to the hard sphere system. We remark that the idea of 
the pair structure being dominated by the hard core^ appears less satisfactory here. 

We now turn our attention to the cluster populations in the Lennard-Jones- WCA systems 
[Fig. EJ^a)]. As before, we consider clusters that are ground states for Lennard-Jones. At this 
higher temperature, compared to Fig. IU(b), relatively little difference is seen between WCA 
and Lennard-Jones, consistent with the concept that in dense liquids, it is the repulsions 
that are responsible for the structure and that attractive interactions have less effect at 
higher temperature. However, the same trend is apparent as was found at the triple point 
[Fig. Hfb)]: Lennard-Jones shows a tendency to form larger clusters than WCA, which seems 
reasonable given its cohesive energy and that these clusters minimise the energy of isolated 
systems. However, as shown by the ratio [N c /N]lj /[N c /N]wca in Fig. [9(a), the difference 
in population is slight. 
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Figure 7: (color online) Pair-correlation functions, (a) Long-ranged potentials: Lennard-Jones 
(red) and WCA (pink) for a well depth of /3elj = 0.4447. (b) Short-ranged potentials: Morse 
(blue) and repulsive Morse (turquoise) according to equation HI Here the well depth /3em = 2.0. 
Cyan denotes the Hard Sphere interaction. 




5A 6A 7A 8B 9B 10B 1 1C 12B 13A HCP FCC 5A 6A 7A 8B 9B 10B 1 1F 12E 13B HCP FCC 

Cluster Type Cluster Type 



Figure 8: (color online) Population of particles in a given ground state cluster. N c is the number 
of particles in a given cluster, N the total number of particles sampled, (a) Lennard-Jones with 
/3slj = 0.440 (red) and corresponding WCA (pink), (b) Morse (po = 25.0) (turquoise) truncated 
Morse (light blue) and hard sphere (dark blue). Note the semi-log scale. 

For the truncated Morse and hard spheres, the cluster populations are not tremendously 
different to the Lennard-Jones case, in fact for smaller clusters the differences are compa- 
rable to those between WCA and Lennard-Jones. In particular, hard spheres show a very 



14 



i— i^i U U U U U U II I I 
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Figure 9: (color online) Ratio of cluster populations in high temperature systems, (a) Lennard- 
Jones and WCA. (b) Truncated Morse and Morse. This plot is the same data as Fig. [8]expressed to 
emphasize the difference between the two systems. Note that in (b) we invert the ratio to consider 
the truncated Morse divided by the attractive Morse potential and plot on a different scale. 

similar population to the truncated Morse, further suggesting that the latter might make a 
reasonable approximation to hard spheres. However, upon adding attractions, the popula- 
tion of clusters drops dramatically. In Fig. Mjo) we plot the ratio [N c / N]tm /[N c /N]m which 
is inverted with respect to Figs. Efa) and |5] in the sense that the truncated system forms 
the numerator and the attractive system forms the denominator. In the truncated Morse 
system, the population of clusters of size m > 10, is at least eight times greater than the 
attractive system. We return to the possible origins of this discrepancy in the next section. 
We note that the attractive Morse g(r) exhibits a split second peak which the truncated 
Morse and hard spheres g(r) do not. Yet, contrary to the notion that the split second peak 
is associated with crystallisation^, in Fig. Wfo) we see precisely the opposite trend: that 
the split second peak is apparently associated with less crystallinity. 

We now plot the population of all identified clusters, noting that it is only for m > 11 
that there is a difference in the ground state clusters for these interaction ranges. That is 
to say, for Lennard- Jones, the ground states are 11C {Ci v ), 12B (Cs v ) and 13A icosahedron 
whereas for the short-ranged Morse the ground states are 11F (C^), 12E (-D3/J and 13B 
(D 5 h). Excluding the attractive Morse system there is rather little variation between the 
different systems. In other words, it appears that for the other higher-temperature systems, 
the topological bond structure may be dominated by the hard core, which is matched in all 
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5A 6A 7A 8A 8B 9A 9B 10A 10B 11A 11B 11C 11E 11F 12A 12B 12D 12E 13A 13B HCPFCC 

Cluster Type 

Figure 10: (color online) Population of particles in a given cluster, at parameters mapped to the 
Morse potential ((3cm = 2.0, po = 25.0). N c is the number of particles in a given cluster, N the 
total number of particles sampled. Here we consider ground state clusters for all ranges of the 
Morse potential 1 --. Note the semi-log scale. 



cases. In particular, although the Lennard- Jones and Morse potentials have different ground 
states, at these high temperatures this is little reflected in the structure. 



IV. DISCUSSION AND CONCLUSIONS 



We have analysed the pair structure and performed a topological cluster classification on 
a range of liquids. The pair structure of Lennard- Jones and longer-ranged liquids is entirely 
consistent with the well-known result that repulsive interactions dominate the local packing 
in dense liquids. Shorter-ranged potentials exhibit a strong response in the g(r) upon the 
addition of attractions. However, one expects that these will be accounted for by the use 
of perturbation theory^. We note that as p — > 0, g(r) — > exp[— (3u(r)] so a short-ranged 
attraction leads to a strong peak at contact as we see in Figure [7(b). 

Although the pair structure of Lennard- Jones and WCA is very similar, we are nonetheless 
able to identify clear differences using the TCC. We find that Lennard- Jones is more able to 
form higher-order clusters than the purely repulsive WCA. These differences become much 
more significant upon cooling to the triple point. Applying the extended law of corresponding 
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states to compare with a weakly interacting longer-ranged system, little effect on the cluster 
population, or pair correlation function is found upon adding attractions. Conversely, in 
short-ranged systems, the radial distribution function is influenced by attractions and the 
cluster population is strongly enhanced upon removing attraction. That we see different 
trends in the short-ranged system is rather curious, and will be investigated further in the 
future. One comment we can make at this stage is that short-ranged attractive systems 
exhibit non-monotonic dynamics as a function of attraction at high densities, in the form of 
a re-entrant glass transition^"— . Whether it is truly appropriate to expect the behaviour 
of short-ranged attractive systems to be similar to long-ranged Lennard- Jones type liquids 
is perhaps an open question. 

We rationalise these three scenarios as follows. The long-ranged Morse is weakly inter- 
acting and compressed. Together, these lead to little response either of the g(r) (the spatial 
distribution of particles) or the topology upon adding attractions. In the Lennard- Jones case 
compression is less important, adding attraction promotes organisation and clustering, how- 
ever the interactions are sufficiently long-ranged that the repulsive core dominates the pair 
structure for both Lennard- Jones and WCA. In the short-ranged case, the hard spheres (and 
presumably the truncated Morse) are close to freezing (here the packing fraction <ft = 0.4875) 
and thus have limited free volume. Adding short-ranged attraction favours configurations 
where the particles are closer to contact, raising the first peak of g(r), and can open up free 
volume. However, considering the second Virial coefficents B\ in table [H these short-ranged 
systems are quite weakly interacting, and there is apparently insufficient cohesive energy to 
promote organisation into clusters. 

Returning to the non-monotonic dynamics of short-ranged systems, one expects that 
the attractive Morse system might exhibit faster dynamics than hard spheres (and perhaps 
the truncated Morse system). Although the hard sphere packing is far from dynamical 
arrest, even so some kind of slowing is expected relative to a dilute fluid. This could 
then be reduced by the short-ranged attraction. Now we have correlated the clusters with 
dynamical arrest 17 and found that arrested states have a high cluster population and that 
it is biased towards higher-order clusters. As a function of density, hard sphere fluids show 
a similar trencU^. Thus we speculate that one possible underlying cause may be related to 
dynamics. However, the hard sphere packing fraction is very close to freezing, and we note 
a substantial quantity of locally crystalline particles in Fig. E^b). Now the colloid-polymer 
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literature 3 ^ would tend to suggest that adding short-ranged attractions widens the fluid- 
crystal coexistence region. Since the hard sphere state point is so close to freezing, and the 
truncated Morse seems similar to hard spheres, it is possible that the equilibrium state for 
Morse (po = 25.0, /fc^/ = 2.0) is crystal-fluid coexistence. It is interesting to note that this 
possibly metastable fluid has a population of HCP and FCC structures around a factor of 
30 less than the stable hard sphere fluid. 

Among the key underlying ideas of clusters in liquids is that they represent energetically 
locally favoured structures^ 2 -^* 3 ^. The most famous of these, the icosahedron, appears only 
in small quantities in this analysis, although it is most prevalent in Lennard- Jones. It would 
be most interesting to investigate whether particles in these clusters are in fact in a low 
energy environment. It would also be interesting to seek a link between structure and dy- 
namics, particularly concerning the recent observation of very different dynamical behaviour 
between the WCA and Lennard- Jones systems 3 ^ and the observation that power-law repul- 
sive interactions seem to recapture the original Lennard- Jones behaviour^. Moreover, other 
mappings have been proposed for example between Lennard- Jones and WCA. Here one can 
place more emphasis upon the dynamics, albeit at some expense in the accuracy with which 
the radial distribution function is matched 41 -. 

Here we have focused on the ground state clusters for each system. Furthermore, liquids 
are by definition at finite temperature, and it may be appropriate to consider the structure 
of clusters at higher temperature in addition to the ground states we have investigated so 
fa r 23 ' — . Conversely, further quenching might favour the ground states beyond the trends 
we have so far seen. Recently, we found we needed around lOfceT of attraction to form 
isolated clusters 44 . 

A final point for discussion is the link between attractions and reciprocal space structure. 
The static structure factor S(q) measures compressibility at wavevector q = 0. While all 
state points sampled show no indication of sny phase transition, one nonetheless expects 
some hint of attractions at low q&£§.. In particular, around the Lennard- Jones triple point, 
for q — > we might expect a factor of two increase in S(q) between WCA and LJ 4 ^ This 
difference in S(q) upon adding attractions is predicted to be most prevalent at a value of 
qo « 1 or around 2ira in real space. This is a rather larger lengthscale than the clusters we 
probe, and in fact Stell and Weis^ show that by qo ~ 3, at the length scales we consider 
here, this effect is much reduced. It would be most interesting to extend the TCC to larger 
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clusters, such that the qa ^ 1 range might be reached. Thus we could directly investigate 
the nature of the change in structure related to this low q behaviour. However, we note that 
the TCC is a particle-based methodology, and may not be sensitive to such delicate changes 
in long-ranged structure. 
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